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THE THEORY OF SOIL CONSOLIDATION 
AND TESTING OF FOUNDATION SOILS 


BY THE DIVISION OF TESTS, BUREAU OF PUBLIC ROADS 


Reported by L. A. PALMER, Associate Chemist, and E. S. BARBER, Junior Highway Engineer 


EVERAL DEVICES used in performing the com- 
pression test described in this report were con- 
structed by the Bureau of Public Roads in 1926 in 
‘accordance with plans furnished by Dr. Charles 
‘Terzaghi. Investigation of testing technique and 
practical utilization of the resulting data have been 
included in the Bureau’s program of research on soils 
for highway purposes continuously since that time. 
As part of this program, the Bureau cooperated with 
the Massachusetts Institute of Technology in research 
on soil mechanics from 1926 until the fall of 1932. Dr. 


Terzaghi was retained by the Bureau as research con- | 


1929. Dr. Arthur Casagrande was employed by the 
Bureau as research assistant and was stationed at the 
\fassachusetts Institute of Technology throughout the 
duration of this cooperative research work. 

Other parts of the Bureau’s general program were 

carried on at the Bureau’s laboratory at Arlington, Va., 
and through cooperation with various State highway 
idepartments and Federal bureaus. 
Special studies have been made at Arlington to 
ijetermine the magnitude of experimental errors result- 
Sug from variables in testing apparatus and procedure, 
nd to determine the agreement between deformations 
wbserved in tests with those computed in accordance 
ith the theory of consolidation on which the tests are 
used, 

he compression test and theory of soil consolidation 
Were described by Terzaghi in 1927.' Two 
Be ports of distinctive value to those concerned with the 
)practical utilization of the compression test data have 
“\so been published.’ 

» On the basis of the Bureau’s investigations of the test 
self, experience in practical use of the test data, and 
eports just referred to, it is now possible (1) to state 
efinitely that the apparatus and testing procedure 
escribed briefly in this report have been found satis- 
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tr at ng that part of the total settlement of soil caused 
pb. @By the loss of water forced vertically out of saturated 
@ompressible soil strata in certain soil profiles. 
vecifically described, or produced by plastic flow, flow 
f water laterally, or causes not discussed, are beyond 
le scope of this presentation. 
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COMPRESSION DEVICE VALUABLE AID IN STUDYING SOIL 
COMPRESSION 


The theory of consolidation of earth materials in- 
olves two problems: (1) Determination of the distri- 
utio. of stresses; and (2) computation of the displace- 
ent. caused by these stresses. In dealing with many 
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ples of Final Soil Classification, Pusiic Roaps, vol. §, ro. 3, May 1927. 

t of the Special Committee on Earths and Foundations, Proceedirgs of 
can Society of Civil Engineers, vol. 59, no. 5, May 1933, and discussion by 
P. Kimball, Proceedings of the American Society of Civil Engineers, vol. 
August 1933. 


Psultant on the work from the fall of 1926 to the fall of | 


other | 


Mactory, and (2) to present a practical method of esti- | 


Settlements of soils in profiles differing from those | 


ering materials and structures the essential | 


problem is to estimate the stresses, since nearly all 
structural materials possess, within limits, the property 
of elasticity. With soils, however, the problem is 
quite complicated for they have both plastic and elastic 
properties; the materials are essentially heterogeneous; 
and water is usually present. If the water is free to 
low from the compressed soil most of the settlement 
results directly from the reduction in voids between 
solid particles, an occurrence that, in this paper, is 
considered apart from elastic and plastic deformations. 

One of the most useful tools in the study of soil con- 
solidation is the laboratory compression device designed 
by Dr. Terzaghi. The compression tests indicate 
within a comparatively short time in the laboratory 
the consolidation characteristics of soils upon which 
depend the rates of settlement and total settlements 
of foundations. Thus, compression tests supplant 
guesswork and “rule of thumb’? methods with reliable 
information for use in the design of foundations. 

Figure 1 iliustrates the essential features of the test- 
ing device used by the Bureau of Public Roads. An 
undisturbed sample of soil with natural structure and 
moisture content, as nearly can be obtained, is 
carefully placed between two porous stones, which 
act as filters during the test. Pressures are then ap- 
plied, ranging from zero to the maximums that are 
expected in service. A burette connected with the 
apparatus as shown in figure 1 affords a means of de- 
termining the permeability of the sample at any time 
during the test. 

The sample during test is representative of a stratum 
of saturated, compressible soil in the earth’s crust, 
sandwiched between two layers of sand or other rela- 
tively permeable material. 

The total amount of consolidation and the speed at 
which it takes place are controlled by: (1) The moisture 
content of the soil; (2) the velocity with which water 
flows vertically through the pores of the sample to the 
filters above and below in the laboratory, or through 
the compressible layer to the sand strata in the field; 
and (3) the frictional resistance of the soil particles to 
consolidation. 

At the instant of load application, the water in the 
compressible sample or layer is considered as furnishing 
all of the resistance to consolidation. Thus, if a suitable 
pressure gauge were connected with the sample in the 
laboratory, it would show that, at the instant of load- 
ing, the water pressure in the sample is practically 
equal to the pressure applied to the sample. It was 
stated during the first conference on soil mechanics and 
foundation engineering, held at Harvard University in 
June, 1936, that experiments in Sweden disclosed that 
the water pressures in compressible understrata were 
approximately equal to the pressures produced by the 
weight of superimposed fill materials just at the time 
construction of the fill was completed. 

As the pressure on the soil continues, water escapes 
from the sample or field stratum, causing the pressure 
to be gradually transferred to the soil particles or 
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values of the pressure, p, greater than about one-tenth 
kilogram per square centimeter, the load-compression 
curve, as shown, is expressed by the equation,® 
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e=1.050—0.442 log p (1) 


\ where ¢ is the voids ratio. The permeability-voids 
\ ratio curve (fig. 2, upper) shows the coefficient of 
permeability of the sample at different stages of 
compression. 
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4 SRARRSRRAARAAR EASES 7/ SR COEFFICIENTS OF PERMEABILITY AND COMPRESSIBILITY 
LLL DETERMINED 


Figure 1.—Device For DETERMINING THE CONSOLIDATION The coefficient of permeability is defined as_ the 
CHARACTERISTICS OF SoILs IN THE LABORATORY. velocity of flow when the hydraulic gradient is unity, 

that is, when the drop in hydrostatic pressure in grams 

. we oe per square centimeter between two points in a soil 
column equals the distance in centimeters between the 

points. The usual practice is to record and plot the 

——I | | | coefficient of permeability in centimeters per second, 

MATED (OASHED LIN as in figure 2. However, for the purpose of simplifying 
a) _— subsequent computations the coe ‘flicient of permea- 
S waa NE bility, as used in this report, is e xpressed 1 in centimeters 
Z 1 : — per minute. The average time-consolidation curve 
“ M (fig. 2, lower) shows the rate of consolidation under 
any given pressure. This curve is the average of the 
N N (APPROXIM QUATION,€ » individual time-consolidation curves for different load 
increments. Essential features of the test procedure, 
the equation for determining the coeflicient of perme- 
ability, and the use of average time-consolidation 

curves have been discussed previously Ky 
ise pe Both the consolidation characteristics of the soil 
| strata beneath a given footing and the vertical pressures 
at various depths produced by the superimposed load 
must be known. In order to solve the many different 
problems involved in computing earth pressures, one 
must be familiar with the mathematical theory appli- 
z AVE OF EXPERIMEN A cable to all types of problems requiring the calculation 
a of earth stresses and their distribution. The treatment 

of special cases is then simplified. 

- ——-s; The coefficient of permeability, ’, is the velocity of 
flow at a hydraulic gradient of unity. If the thickness 
of the soil sample, with a voids ratio of e, is d, and if the 
voids ratio could be reduced to zero, the thickness 


Figure 2.—PERMEABILITY, COMPRESSION, AND Conse.ipation | would become dp, that is, 
CURVES FOR A SAMPLE OF COMPRESSIBLE SOIL. 
d=d,(1+e)-_-- 2 
“soil skeleton.”’” Eventually, equilibrium is reached 
and further consolidation ceases when the water pres- 
sure becomes equal to zero and the soil skeleton pro- 
vides all of the resistance to consolidation. 

The sample in the compression device is subjected to 
pressure that is increased by increments. Each in- 
crement of pressure acts until apparent equilibrium is 
attained before the next increment is added. Water kd “3 
escapes under pressure, entering the porous stones and | ky dy 
flowing from the outlets as indicated in figure 1. that is, 

Consolidation characteristics and permeability of the 
laboratory specimen are recorded and the data are ky cle a | 
plotted as shown in figure 2. The load-compression | 
curve shows the voids ratios at which the soil sample | , oe ciple sof Final Soil Classification, by Dr. Charles Terzaghi, Pupuic Roa 
attained apparent equilibrium under different pressures. | ‘When p dotted tom seanltogualtenia seabe, 0 helen uiitned to the natuesl exis. 
The voids ratio is defined as the ratio (by volume) ee ce gape A <p tego pov sag the an ies 
of voids to solids for a given sample of soil. curve is 

Apparent equilibrium in the sample is reached in a | e=B-Z log p 
relatively short time, and for all practical purposes is | Where B is the intercept of the straight-line portion of the curve at the 1.0-kilog? 


me 7 2 per-square-centimeter ordinate and — Z is the sl »f the log plot of tt 
true equilibrium. Actually, however, the attainment | range in value ofefor unit changeoflogp.. 


‘A Method of Predicting Settlement of Fills Placed on Muck Beds, by | 


The rate at which water would permeate the soil 
sample if the difference between the em water levels 
were do, rather than d, is designated as hp. It has been 
shown ' that the rate of permeation increases as the 
hydraulic gradient increases so that from this fact and 


equation (2 
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of true equilibrium would require infinite time. For ! oveson Pesuc Rorvs, vol. 16. he 12, Pebeuney tose 
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It is convenient to take e as the average voids ratio of 
a soil sample compressed from an initial voids ratio, 
¢;, to an equilibrium value e:, so that 


€; + €2 
») . aa 
On substitution of this value for e in equation (4), 
one obtains 


k, (6 


The coefficient of compressibility, a, is the average 
decrease in voids ratio per unit increase of pressure. 
If the initial voids ratio, e;, is the condition at equilibrium 
under a pressure of p, grams per square centimeter and 
¢, is the corresponding condition at a pressure p, grams 
per square centimeter, then 


The values of both 


ratio, but at different 


a and & decrease with the voids 
rates. However, the change in 


the value of the ratio, ; The ratio, 
( 


>is fairly small. 


/ 


is called the coefficient of consolidation, ¢. From 
equations (6) and (7) there results from substitution for 
al dl , 
} p ) 2k 
( Ss / \% 
a ( ( 2 ¢ r 


The itility of this equation Mia be illustrated by use 
of the data plotted in the upper curves of figure 2. 

From this curve, the pressures, P; and Pe, correspond- 
Ing to ¢ 1.00 and e.=0.80, are 1,300 and 3,680 grams 
per square centimeter, respectively. The average voids 
ratio is 0.90 and from the permeability curve, the 

efficient of permeability corresponding to a voids 
tio of 0.90 is 4.6% 1078 centimeters per second 


; 3.680 1.300 92x 19°? re 17 
Then ’ ~ 60—0.0173 
1.00—0.80 “*2+1.00+0.80’*> as 
Ci Ineter per minute. 
“Similar computations can be made, taking different 
values for é, and e, and the corresponding pressures, 


] nd p. from the curve, finding the coefficient of 
permeability corresponding to the average voids ratio 
as lllustrated above, and substituting these values in 
equation (S By this procedure, the values of ¢ given 
le 1 are computed. These values are plotted 
the corresponding voids ratios in figure 3. 


in t le 


against 


for voids ratios of 0.6, 1.2, and 1.4 were computed from 


the relation, 
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This paper introduces only the most essential steps 
in some of the mathematical derivations used both for 
the computation of earth pressures and for computation 
of formulas applicable to the consolidation of relatively 
Impervious, saturated strata caused by these pressures. 
lor readers more concerned with direct application of 
the fundamental principles and formulas than with 
their development, a discussion of the use of the theory 
is presented in detail in part 1 of the paper. The mathe- 
matical derivations are included as part 2 for those 
interested also in the theoretical development. 


!! 


The manner of applying the test data in practice | 
varies, depending upon the pressure distribution condi- | 
tions beneath the loaded area, the thickness of the | 
compressible layer, and whether it is free to drain both | 
at top and bottom, as in the compression test, or from | 
only one of the bounding surfaces. 

When a load is applied on the top of a compressible | 
soil stratum, conditions may be such that the pressure 
resulting from the load is the same from top to bottom 
of the stratum, that is, the pressure distribution is 
uniform. On the other hand, the conditions may be 
such that the pressure varies, i. e., increases or decreases 
as the depth increases. In this paper the pressure 
change per unit of depth is considered to be either zero 
or constant. If it is zero, the pressure distribution is | 
uniform; if it is constant, the pressure distribution is 
either triangular or trapezoidal. The cases of pressure 
distribution are illustrated in figure 4 and are as follows: 

Case 1.—Rectangular pressure distribution. The 
pressure is uniform from top to bottom of the compres- 
sible stratum. 

Case 2.—Triangular pressure distribution. The pres- 
sure varies from zero at the upper bounding surface of 
the compressible stratum to a finite value greater than 
zero at the lower bounding surface. 

Case 3.—Triangular pressure distribution. The pres- 
sure varies from a finite value greater than zero at the 
top of the compressible layer to a pressure of zero at the 
bottom. 

Case 4.—Trapezoidal pressure distribution. The 
pressure increases with increasing depth in the com- 
pressible stratum. The pressures at both top and 
bottom surfaces are finite values, each greater than zero. 

Case 5.—Trapezoidal pressure distribution. The 
pressure diminishes with increasing depth in the com- 

ressible stratum. The pressures at both top and 

ottom surfaces are finite values, each greater than 
zero. 


FORMULAS DEVELOPED FOR VARIOUS TYPES OF LOADED AREAS 


The pressure distribution throughout the compres- 
sible layer depends on the shape and magnitude of the 
loaded area at the ground surface and on the load 
distribution over this area. In the examples cited in 
this paper the load distribution is taken as uniform 
over the loaded area. 

The five types of loaded areas considered are: 


A.—A loaded area of great length and width. 

B.—A spot load. 

C.—A relatively short and narrow, rectangular 
loaded area. 

D.—A relatively narrow loaded area of great 
length. 

E.—A circular loaded area. 


The pressure per unit area at any given depth below 
a loaded area of great extent is taken as equal to the 
unit pressure on the loaded surface plus the weight of 
superimposed earth above the point in question. No 
special formula is required for computations in this 
simple case. 
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PART 1.—APPLICATION OF THEORY OF SOIL CONSOLIDATION 
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Figure 4.—Five Cases or Pressure DISTRIBUTION IN Sort 
STRATA. 
where 3 
» 
» o 
K- = r vy ieee 10 
2x(1 } =) 
P=the spot load; and 
p,—the verti al pressure per unit 


area at a vertical distance, 2, 
below a spot load on the sur- J} 


j 


face and at a horizontal dis- @ 

tance, r, from the load. . 

For any distance on the vertical axis below the load, 

9 Se 

° = » 4 

ris zero and A=; : 
Equations (9) and (10) can be used to determine the & wl 


vertical pressure produced by any load, uniform) 

distributed over a rectangular area, provided the ratio 
of the longer side of the rectangle to z, the vertical § 
depth to the point at which p, is to be computed, is @ 
not greater than one-half. This condition is illustrated bs , 
in figure 5. If, however, this ratio exceeds one-hal!, § 


any 






then the procedure is to subdivide the rectangle into ® Jin, 
rectangular elements (all of equal area) such that the® 
| ratio of the longer side of a single element to the ver-§ 
tical distance to the point in question is not greater wh 
than one-half. If the point at which it is desired to, 
know the vertical pressure is at a very considerable & = 
depth below the loaded footing, then the procedure o! 
subdivision may be unnecessary even in the case of§ 
large loaded «areas. For example, the base of thei 
Washington Monument is 125 feet square. At any per; 
DF t 
point 250 feet beneath the base the ratio is ea 18 Sh 
a - sibli 
and hence the entire loaded area may be considered a3 wat, 
a point load at this depth. me Th 
When subdivision is necessary it is considered that 9 filte: 
the load on each element is concentrated at its center. §§ pres 
| To determine the unit pressure on the vertical center Tedu 
axis, the horizontal distance r from the center of each] lon 
| element to the center of the footing is computed and ee 
| ‘ { 
the ratio, os is evaluated. This is substituted in the ie Say i 
Z me Can \ 
expression for K, equation (10), and K is computed fo the 


To compute the pressure below a spot load, use is 


made of the equation, 


"iE 


Having done this with all of the 9% dimi, 


the element of area. 


different elements of area, the K values are addei ji the ¢ 

together and the sum is substituted in equation (9), th = Th 

value for P in this equation being the total load on eaci fy Comp 

individual subdivision of area. This procedure will b’§ oad~ 

further illustrated later. For s 
For a relatively narrow loaded area of great lengtlf 

take a vertical cross section of the loaded strip. \t') 

point on the vertical center axis, p, is computed fro! 

the equation, for so 
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hiGURE 5.—RECTANGULAR AREA LOADED UNIFORMLY. (AREA 
or Type C, ILLUSTRATING THE PRINCIPLES OF Botu TyPpEs 
B anp C., Z 
~B 
ae 
es og | are (11) 
“ 4 T 
4 where a= twice the angle, in radians, whose tangent 
is half the width of the strip divided by 

the vertical distance z to the point in . — , . m ¥ 
\ ; 4: 7 . ar a . ‘ ‘IGURE 0.--LYPES OF sOADED AREAS. PPER, ERTICAL 
ul 3 question (see fig. 6); and p= the load Cross-SECTION OF A STRIP OF GREAT LENGTH (LOADED 
si per unit area of the footing. ArEA OF Type D). Lower, a Crrcutar LoapED AREA 
dg To compute the maximum vertical pressure, p,, at cess 5). 

, any point below a circular footing and on the center 5 
+ spe fi ) : > oti ; l €2 
oF line (see fig. 6), use the equation, Q | ¢ “|p, ah te Se et oe 1 

e * \ i 
~ | 1, = p (1—cos*B)__- aad _._(12) . 

i P=! i P where e;= voids ratio of the sample, taken as the aver- 
Ti where Pp the load per unit area of the footing age volds ratio In the compressible 

} 5 eo > tr 10Pr hd > : 
and 8=the angle whose tangent is — (see fig. 6). stratum prior to the application of addi- 
: Z . tional load at the surface; 

ATE OF SETTLEMENT COMPUTED BY SIMPLE FORMULAS ¢2= voids ratio of the sample taken as the average 
of voids ratio in the compressible stratum 
he] ii the compressible layer is sandwiched between two upon reaching equilibrium after the appli- 
ny §} permeable sand-gravel layers, one above and one below, cation of additional load at the surface; 

1 J} water can escape in two directions. If the compres- d,=thickness of soil sample at apparent equilib- 

2) sible layer is overlaid with sand and rests on rock, rium at voids ratio, e;; and 
as} Water can escape from but one face of the stratum. D,=the thickness of the soil stratum at voids 

© There is a double filter in the first case and a single ratio, é;. 





sat (filter in the other. In either case, under a given 
er. pressure the thickness of the compressible layer is 
ter fm Teduced to the same extent, ultimately. However, the 
ach @ Tates of reduction in thickness of the compressible 
und layer will obviously be different in the two cases. Con- 

> sequently, the settlement during a relatively short time, 
the Say i. 6 months or a year, will be greater if filtration 
_ § can o-cur at both the upper and lower boundaries of 
lo!) the compressible stratum. The filtration effects a 
the diminiition in thickness and the earth’s surface above 
dei the compressible stratum settles. 
th) The total amount of vertical consolidation Q is 
acife COMputed by substitution of the data shown by the 
| be load-c: npression curve in the following expressions: 

© For soil sample, 
oth fi 
\ Q -[ oe Ja eo (13) 
role ney 

for soi stratum, 








The rates at which soil strata settle are computed 
from test data and field observations by simple formulas. 

The theoretical time-settlement relation for soil 
samples or strata sandwiched between two porous 
layers and with either rectangular, triangular, or 
trapezoidal pressure distribution or for strata with 
_ rectangular pressure distribution and with a porous 
layer only at the upper face is given by the expression, 


ke; . (15) 


where 9g,;=percentage of ultimate consolidation at the 
end of a definite period of time, t, of load 

application; 

e= Naperian base. 

The exponent, N, is related to the thickness 
of the soil sample, the period of load appli- 
cation, t, and the speed of egress of water 
| from the soil sample or stratum. 
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Experience bas shown that the time-consolidation 
curves furnished by the compression test on compressi- 
ble soils (see fig. 2, lower curve) conform in shape with 
the theoretical curve of equation (15) from g, equals 
zero to g, equals 90 percent on the theoretical curve. 
The experimental values of q, attributable to per- 
meability alone, are usually from 60 to 100 percent of 
the actual ultimate consolidation as obtained in the 
laboratory with the given loading conditions. Only 
laboratory values of g; within the range in which the 
experimental curve agrees with the curve showing the 
theoretical relationship are valid for computations 
involving the use of equation (15). Beyond the limit 
of consolidation as effected by the egress of water, 
further consolidation of most soils is caused by particle 
rearrangement and does not result from permeability 
as assumed in the derivation of equation (15). 


TIME-CONSOLIDATION RELATIONS GIVEN FOR UNIFORM, TRIANGU- 
LAR, AND TRAPEZOIDAL PRESSURE DISTRIBUTIONS 


Considering a laboratory sample, the equation for 
the exponent N is 


When a definite numerical value is assigned to q,, NV 
has a numerical value, designated as \,. The expres- 
sion then becomes 


N( %) 


Baie (16 
2.47¢ 


where =period of load application, in minutes. 
dy=thickness in centimeters that the labora- 
tory sample would have if it were possi- 
ble for the voids ratio to equal zero. 
d,—twice the maximum distance water travels 
when flowing from a sample of dy) thick- 
ness sandwiched between two filters. 
c=coefficient of consolidation of compressible 
soil, in centimeters per minute. 
To express t of equation (16) in years instead of min- 
utes, and d, in feet instead of centimeters, the following 
computations are necessary: 


1 centimeter ne (3 feet 
and 
1 minute — —_ vears 
365.24 * 24 « 60° ‘ 
Then 


w(4) 


N. $) 
ten NSS (2.5412)? | 
~ 2.47¢ ~365.24K24X60 1,400 ¢ 9°": 
In this equation d, is in feet. 

For a compressible stratum in the field, sandwiched 
between two permeable layers of soil, the time required 
for a given percentage of consolidation to take place is 
computed from the equation, 


1,400 ¢ ° 


| 
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where /),—the thickness in feet of the stratum at ¢—0. 

When a compressible soil stratum rests on rock and is 
covered with a porous layer, and when the pressure dis- 
tribution is uniform, the time-consolidation relation is 


N,D,? 1g 
== vears ) 
1,400 c* 

When a soil stratum rests on impermeable rock and 
is covered with a porous layer, and the pressure dis- 
tribution is triangular (case 2), the expression becomes 


ND? 


(19) 
1,400 ¢ ‘ 


years 


where V, is a definite numerical value corresponding 
to an assigned value for q,, the percentage consolida- 
tion at the time, f, for this tvpe of pressure distribution 
and boundary conditions. For values of g, and q., the 


corresponding numerical values of NV; and \V, are given 
in table 2. 
For triangular pressure distribution (case 3) with a 


permeable layer only at the upper bounding surface of 
the compressible stratum, the time-settlement relation is 


N;D,° 
1,400 ¢ 


years ee! 20 


For values of g, the corresponding numerical values o! 
N, are given in table 2. Here 3 denotes the percentage 
of consolidation at any time, ft, with this type of pres- 
sure distribution (case 3) and with the stated boundary 
conditions. 

The procedures followed in computing the theoretic:! 
values of N,, V2, and N3, corresponding to values of 








M1, Jz, and gq; will be shown on pages 7 and 8 preceding 
the illustrative examples. 
Ped 
ot 
von 
T ABLE 2 J alues of N . N . and N; correspond ng to value oa 
’ ! 2 / & 
for cases 4 3, respective ny 
91, G2, and qz for cases 1, 2, and 3, respec i 
Numerical values of 
Consolida 
tion qi, 92, 
i N; \ \ 
\ ce at 
U0 0.00 0.00 OU 
OK OO) . WOR ; 
02 12 On bs 
] O4 Is 01 | AES 
20 OS 25 i ke 
25 12 a 03 Fe 
0 17 v On . 
35 24 47 08 Be 
40 3 7] 12 by 
45 . 39 63 et 8 
50 19 73 24 c4 
5 59 84 .32 ; 
6 aa 95 42 : 
6 84 1, 08 54 a 
70 1, 00 1. 24 . 69 | 
75 1, 18 1. 42 8S Fim 
SU 1. 40 1. 64 1, OS 8 
85 1. 69 1. 93 1. 36 i 
90 2.09 2. 35 1.77 Lae, 
95 2. 80 3. 02 2. 54 2 
100 Infinity Infinity Infinity m4 
ey 


With trapezoidal pressure distribution, the pressw] 


resulting from the applied load may either increas 





(case 4) or decrease (case 5) with depth, and has @ 
finite value at both the upper and lower boundaries (| 


surfaces of the compressible stratum. The ratio o! thi 


pressure at the upper boundary to that at the lower} 
The value of w is unity for union 


designated as w. 


| and either zero or infinity for triangular pressure dj 
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tribution. For trapezoidal pressure distribution, (case 
4), wis less than | and greater than zero. An interpo- 
lation factor, J, is applied in this case in computing the 
rate of settlement if the compressible stratum is overlaid 
with sand and rests on rock. When this condition 
obtains, interpolation is between the two extremes of 
uniform and triangular pressure distributions. With 
the compressible stratum bounded below by sand and 
above by rock, and with a trapezoidal pressure distribu- 
tion such that u exceeds unity, (case 5) an interpolation 
factor, J’, is used in computing the rate of consolidation. 

The expressions for J and ./’ in terms of w are as 


follows: 
Ee 2)u- ‘] 
log 
J ' w( |] +-w) 


log. 

| 99 
log z 
2\n—2Z) 


u(r—2)-4 
log} “\7— 
. (7 2iie-7 
J 
Values for J and J’ corresponding to values assigned to 


vw are listed in table 3. 


Kor case 4, . 


— IC 


For ease a , 


The time-settlement relations for cases 4 and 5, 
trapezoidal pressure distributions, are: 

N,D,? 
( use 4 "ears 28 

iS¢ . 1.400 ¢ Veal “ 
. . N;D,? ; 

use o rears - 

: 1400¢ %°™ -“ 

where N,=N,4+- dJ(N2—N, (23 
ind N;=N3+J’(N,—N) - 26 


L ABLE 3. Values of u, J and J’ 
ttlement for the twot yy pes 


ised in computing the rates of 


of tra pe zoidal pressure distribution 


Pressure increasing with Pressure diminishing with 
depth, u varying from 0 depth, u varying from 1 
to 1 to infinity 
Value of u Value of J Value of u Value of J 

oo 1 0OO 1.0 OO 
12 92 

l st 1.4 Mf 
1.6 so 

4 Hu LS iF 
1.Y ‘ 

r 2. 0 71 

a & 67 

{ if} 2. 5 H2 
$0 55 

3.5 w 

1.0 4 

h 2 15 42 
5.0 iy 

7 ; 6.0 ; 
7.0 ) 

S 12 s.0 27 
9.0 2 

} Of 10.0 z 
12.0 1) 

1.0 00 15.0 17 
20.0 13 

100.0 02 


FAIRLY CLOSE AGREEMENT FOUND BETWEEN THEORETICAL AND 


TEST VALUES 

lhe time-consolidation curve (fig. 2) represents the 
tverugve curve for pressures of p=1.079, 2.118, 4.196, 
and S.352 kilograms per square centimeter. The time- 
Consolidation curve for the load p=0.61 kilogram per 
square centimeter was obviously influenced by experi- 
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mental error and therefore was not included in com- 
puting the values for the average curve. 

The divergence of the average curve from the curve 
given by equation (15) may be determined by means 
of equation (16) as follows: 

For different values of g,;, say 25, 50, and 75 percent, 
the values of ¢ are obtained from the curves. The 
average values of e at the different stages of loading 
are also obtained from these curves. Then the values 
of c, corresponding to the average values of ¢, are read 
from the curve in figure 3; the value of d, is computed; 
the numerical values are substituted in equation (16); 
and the values for NV, may thus be computed, assuming 
that the relationship in equation (16) obtains. 

The thickness of the sample, d, at a voids ratio, 
¢=1.670, was found in the test to be 1.130 centimeters. 
Therefore the thickness d, when e=0 would be 0.423 
centimeter. 

The values for .\V, thus obtained (see table 4) for 
equal values of gq, are averaged and found to be 0.11 
for q,=25, 0.44 for g,=50, and 1.30 for g,=75. 

The theoretical values of .V, for the same values of 
q, are shown in the last column of table 4. The agree- 
ment between the determined and theoretical values is 
fuirly good. The procedure in computing these theo- 


retical values of .V, follows: 





| BLE 4 Val es Jor q &,c, and \ obtained f ‘om laborato y data 
\ ue f Value f 
( ; iW} Nit 
¢ g puted 
rage fr ita 
A I 

O10 t é 

O79 Ollft | 44 

» O140 47 » 

( 1 2 

. O159 42 44 

4 O1l6t 20 s 

x7 { 74 12 2 

2 x 4 s 180 42 at) 

% O83 2 s 

CSS { 

8 O1R9 14 14 

} 22 Is 

' 
For 25 percent, average .\ mputed from data in column 5)=0.11. For 
0 percent, average N ited from data in column 5)=0.44. For 75 per 
t, average .N; (computed from data in column 5 1.30. 


COMPUTATIONS GIVEN 
TRIBUTION 


FOR VARIOUS TYPES OF PRESSURE DIS- 
AND BOUNDARY CONDITIONS 


Case 1.—The computations for case 1 cover uniform 
pressure distribution, irrespective of boundary condi- 
tions, and triangular and trapezoidal pressure distri- 
bution when the compressible stratum is bounded on 
both top and bottom by permeable strata. To compute 
the value of \, corresponding to any definite value as- 
signed to g; in the expression, 


(15) 


one may, as an approximation, 


except the first. Then 
8 
qi=1——<e (27) 
or 
8 
l—q=—e N . : (28) 


Ee 


" 1 . 
For g,=3 equation (28) becomes 


) 
Ss : 
Li ri) ~N 1 
z ¢ 


and 


1 S : 
log ¢ 5 log ¢ =2—M, 


where “log ,”’ designates the natural logarithm. 
s 1 8 1 
N,=log . =—log , : -2.3026( log = —log ) 
l DS é r- ? € 9 DS Tr S 2 


where ‘‘log’’ denotes the base 10 system of logarithms. 
Hence .V,=0.48. 

The more accurate procedure would be to assign 
values to .V, include several terms of the series, and 
solve forg,;. In this manner it is found that for V,=0.49 
the value of q, is 0.50 or 50 percent. (See table 2.) 
However, it is seen that by neglecting all terms of e 
except the first in equation (15), the error, 0.49—0.48 
=(.01, is small. 

Case 2.—Triangular pressure distribution with zero 


| 
| 
| 


| 
| 
| 


pressure at the upper boundary, when the compressible | 


stratum is bounded by an impermeable stratum below 
and a permeable stratum above. For this condition 
the theoretical time-settlement relationship is given by 
the expression, 


] = 

q=1—-;(e *——e "+e 7" —, + ...)--.(29 
1 al 24 +t 125 29) 
By neglecting all terms of e except the first, one may 
write, 


32 
1— > : - oi 30) 
2 = (3 
l . 32 ao 
for Q=s3’ N, loge 3—loge5 =Q.72. 


This compares favorably with the value, N,=0.73 for 
q2=50 percent, table 2, obtained by taking more terms 
in the series than the first in equation (29). 

Case 3.—Triangular pressure distribution with zero 


pressure at the lower boundary, when the compressible | 


stratum is bounded by an impermeable stratum below 
and a permeable stratum above. It will be shown in 
part 2 of this paper that 
q3=2q1— 42- - _(31) 
The simplest procedure in computing values of N, that 
correspond to values of g; is to take the case, N;=N, 
=N;. For Ni =N,=—N3;=0.12, it is seen in table 2 that 
1=0.25 and q,.=0.10, so that g,=20.25—0.10=0.40, 
hat is for g,=40 percent, N;=0.12. 

Case 4.—Trapezoidal pressure distribution with pres- 
sure increasing with depth, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. The value, N,, is 
computed from equation (25), which is Ny=N,+ 
J (N.—N,), applicable for q=q,=q. J in turn is 
computed from equation (21) for any given value of w. 
For example, when u=0.50, J=0.36 (see table 3). For 
n‘h=g=Uu=50 percent, VN, =0.49 and N,=—0.73. Then 


N;=0.49 +0.36 (0.73 —0.49) =0.58. 
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When this value is substituted for N, in equation (23), 


N.D¢ peg 
"ears 
1,400 ¢- : 
the time required to produce 50 percent of the ultimate 
The value of ¢ is 
obtained from laboratory data according to the proced- 
D 
= "a as already shown 
+e 
(equation 2). The thickness, D,, of the compressible 
stratum prior to applying a load at the surface is known 
from the test borings. If ),=20 feet and e;=1.5, then 


ure already indicated and Do= 


20 : - , ‘ 
Dy=5-==8 feet. And if c=0.01 and N,=0.58, then 
2.0 
0.58 & 64 - ee , : 
t=— id “= =2.65 years. This is the time required for 


50 percent of the ultimate consolidation to take place 
under the given conditions, since a value of Ny=0.58 
corresponds to 50 percent consolidation. 

Case 5.—Trapezoidal pressure distribution with pres- 
sure decreasing with depth, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. For q,;—4q;=4q;, the 
value, V; is computed from equation (26), which is, 


N;=N3+J’(Ni—Ns), 


for any given value of vu. From table 3, J’ 
0.49 and N;= 


0.71 when 


u=2. For q.=9q:=4q;=50 percent, V,= 
0.24 from table 2. Then by substitution, V;=0.24-+ 
0.71(0.49—0.24)=0.42. Now take the conditions al- 


ready assumed for case 4, that is, c=0.01, e;=1.5, the 
average voids ratio in the stratum prior to loading, 


and D=20 feet. Then 
Nf’ ., 0-42 x64 - 1.9 vears 
1.4000." 14 or 1.3 years. 


DISCUSSION OF CONDITIONS GOVERNING ILLUSTRATIVE EXAMPLES 
THAT FOLLOW 


In the following examples it is assumed that th: 
material composing the compressible stratum has the 
consolidation characteristics derived from computations 
of the data plotted in figures 2 and 3, and that the 
weight per cubic foot of this material as well as that 
of the sand above it is 100 pounds. Cases 1 to 5 
inclusive are linear types of pressure distribution, that 
is, the pressure distribution diagrams are composed of 
straight lines. All of these five types of pressure <is- 
tribution give the same time-consolidation relation, 
equation (17), when the compressible stratum is 
bounded both above and below by drainage cours:s 
(boundary condition A). Equation (18) is the time- 
consolidation relation when there is but one drainage 
course (boundary condition B) and the pressure distt!- 
bution is rectangular (case 1). When there is but one 
drainage course, the stratum resting on rock and beiig 
overlaid with sand, equations (19), (20), (23), and (24 
are the time-consolidation relations for cases 2, 3, 4, 
and 5, respectively, and computations are made wiih 
the aid of tables 2 and 3 and the interpolation formu!:s, 
making use of J and J’, that is, equations (21), (2.), 
(25), and (26). 

Cases 1 to 5 inclusive may be defined in terms of /, 
the ratio of the pressure at the upper surface of 1 .1¢ 
compressible stratum to the pressure at the lower sur- 
face, that is: 
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Case 1, rectangular pressure distribution, w= 1. 

Case 2, triangular pressure distribution, w—0. 

Case 3, triangular pressure distribution, «—infinity. 

Case 4, trapezoidal pressure distribution, ~ varies 
from 0 to 1, 

Case 5, trapezoidal pressure distribution, uw varies 
from 1 to infinity. 

In solving any problem involving case 4, the formulas 
and principles of cases 1 and 2 are applied. When the 
pressure distribution is of the type case 5, application 
is made of the formulas and principles that pertain to 
eases | and 8. Cases 1 and 2 are both limiting condi- 
tions of case 4, and cases 1 and 3 are both limiting 
conditions of case 5. The types of pressure distribu- 
tion most frequently encountered in foundation prob- 
lems are cases | and 5, more especially 5. 

Example 1.Case 1 pressure distribution; type A 
loading. Consider a uniform load of 3 tons per square 
foot over an area of great extent. Immediately be- 
neath the loaded area there is a laver of sand, 10 feet 
thick. Below the sand is a compressible stratum also 
10 feet thick and below this stratum there may be 
either (A) sand or (B) rock. It is desired to know the 
total settlement and the time required for 25, 50, and 
75 percent of the total settlement to take place. 

Use equation (17) for boundary condition A and 
equation (18) for boundary condition B for the time- 
consolidation part of the problem. It is assumed that 
consolidation of the compressible stratum resulting 
from its own weight and that of the sand above it has 
reached completion prior to application of the surface 
load. 

Prior to loading, the pressure is 1,000 pounds pet 
square foot at the top and 2,000 pounds per square foot 
at the bottom of the compressible layer. In figure 2, 
these pressures correspond to voids ratios of 1.28 and 
1.06 respectively. A pressure of 1 kilogram per square 
centimeter is very nearly 1 ton per square foot so that 
the units of pressure of figure 2 (load-compression curve 
are directly applicable. The average voids ratic prior 
to loading is then 

1.28 


”) 


1.08 : 
mi €}. 


\fter loading, the pressure at the top of the stratum is 
3.5 tons and at the bottom, 4 tons per square foot. 
These pressures correspond to the voids ratios (see fig. 
2), 0.80 and 0.78, respectively, the average, ¢. being 
0.79. The total settlement is computed by equation 


14), that is, 
é,—< 1.17—0.79 : 
0 pee if _— O=1.75 feet. 
oie = Thai *! le 


there were no voids whatsoever in the material 
composing the compressible stratum, its’ thickness 
Wi ld be 


l 
dD), =X 10 


l1+e, 2.14 


D), 4.61 feet. 


The average voids ratio is reduced by the surface load- 
ing from e;=1.17 to @=0.79. Within this interval 
1.17+-0.79 


“0.98. 


there is an average e of From figure 3, 


& voids ratio of 0.98 corresponds to a value for c, the 
cocifieient of consolidation, equal to 0.0157. 
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Now apply equation (17) for boundary condition A, 
the compressible stratum sandwiched between two 
drainage courses. The times required for 25, 50, and 
75 percent of ultimate consolidation (1.75 feet) to take 
place are, using table 1 for values of Nj, 


vi(3) 


Ls 
1.400 ¢ 


fol q 25 percent, and \ 


0.12) ( st 


0.03 year or 1] days, 


1.400 0.0157 
0.49 a ) 
for g,=50 percent, .V,=0.49, f 1.400 » 0.0157 
0.12 year=44 days, 
1.18 bo) 
for q,=75 percent, N;=1.18, f 1.400 » 0.0157 


0.285 vear=104 days. 


lor condition B, the stratum covered with sand and 
resting on rock, the values of N, corresponding to 
values of gq; are the same as for condition A. Solving 


3 18) f u tf 0.12 « (4.61)° 0.12 
eque nh | )tor g 20 ‘reent, ? —o AZ 
jula on OI ] pel en 1.400 20.0157 

vear=44 days. Similarly, for q,=50 percent, t=0.48 


year, and for q,=75 percent, t=1.14 years. 

Twenty-five percent of 1.75 feet, the ultimate settle- 
ment, is 5.2 inches. For condition A, this amount of 
settlement occurs in 11 days; for condition B, 44 days 
are required for this amount of settlement. 

Example 2.— Case 5 pressure distribution; type C 
loading. Instead of a loading area of great extent as 
in example 1, there is a rectangular rigid footing, 10 
feet wide and 30 feet long, carrying a uniform load of 
3 tons per square foot. ‘This is at the ground surface 
above the same compressible stratum described in 
example 1. The average voids ratio, ¢,, of the com- 
pressible layer prior to loading is taken as 1.17; its 
thickness prior to loading is 10 feet; and the weight of 
the sand at the upper level of the stratum is 1,000 
pounds per square foot 

The vertical distance, z, of the footing above the 
upper boundary of the stratum is 10 feet. The ratio 
of the length of the footing to zis 3. This ratio is too 
large to warrant taking the uniform load on the footing 
as a point load at its center. For purposes of computa- 
tions it is necessary to consider the footing as sub- 
divided into rectangles such that the ratio of the length 
of a single rectangle to 2 does not exceed one-half (see 
figs. 5 and 7). The footing may be divided into twelve 
5-foot squares to meet this requirement. 

The total load on the footing is 3 1030900 tons. 
The load on a single square is 75 tons and it is con- 
sidered that this load is concentrated at the center of 
the square. The horizontal distance from the center 
of a square to the center of the 10- by 30-foot!foundation 
is r and this is found from the properties of a right 
triangle. For example, the center of square number 5 
from C, the center of the footing, is the square root of 


5\° 5\? 5 Tr ‘ ‘ 
( >) +(2) ‘or ~==r. The same value for r is obtained 
2 2 /2 
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Figure 7.—ReEcTANGULAR LOADED AREA (TyPpE C 


. = rr : r . ; 
for squares 6, 7, and 8. The ratio, —+ appearing in 


. . " 
equation (10) 1 or 
"102. 2y2 
These values and those obts cilia for the other squares 
(fig. 7) are given in table 5, together with the value of 


; for squares 5, 6, 7, and 8 


, . . : . % 
K obtained by substituting the value for — in equa- 


tion (10), 


9 
o 


The three values for K, corresponding to the three 
groups of load numbers given in the first column of 
table 5, are added together, their sum being 0.5404. 
This must be multiplied by 4 to give K for the entire 
footing. The vertical stress at a point 10 feet below 
C and on the vertical center axis through C is then 
calculated from equation (9). That is, 


; KP 75 
P:~ 2 ~ 100 


K 


x4 xX 0.5404 = 1.62 


? tons per square foot 


TABLE 5. 


-—— 


Computations of vertical stresses below a rectangular 
footing 10 by 30 feel 


7 ? Coeffici 
Load numbers (see fig. 7 Distance, 7 Ratic of "Kk sie 
| 
| Feet 
. 5 
| 5, 6,7, and § é = =(), 354 0. 3555 
v2 ev 2 
‘ 0 10 
3, 4, 9, and 10 x. : ; 0.790 1420 
¢ 26 
1, 2, 11, and 12 v y ; O74 0429 
Total 5404 


By taking z=20 feet and computing in the same 
manner, it is found that at the lower boundary of the 
compressible stratum and on the vertical center axis of 
the footing, p, is 0.73 ton per square foot. The pres- 
sure is a maximum for points along the vertical center 
axis of the footing and the distribution of this maximum 
pressure is trapezoidal, case 5. At points removed 
from the vertical center axis there is also case 5 pres- 
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FIGURE 8 DISTRIBUTION OF MAXIMUM PRESSURE UNDER 
RECTANGULAR Footrina. Sor LAYER CONTAINED BETWEEN 
Two PERMEABLE SAND LAYERS. 


sure distribution, but the value of p, at any level 
less than the value of p, for the same level on the 
vertical center axis For computations consider only 
the maximum pressures, which are along the vertical 
center axis 

The shaded trapezoidal area marked “inactive” in 
figure 8, denotes the pressure distribution resulting from 
the weight of earth material alone. The trapezoidal 
figure marked ‘‘active’’ denotes the distribution of max- 
imum pressure in the compressible stratum that results 
from the uniform load on the footing. It is again 
assumed that consolidation of the stratum caused by 
its own weight and that of the sand is complete at the 
time the surface load is applied and hence all of this 
pressure has been transferred from water to soil and 
therefore cannot act in combination with the applied 
surface load in causing consolidation. The applied load 
producing the active trapezoidal pressure distributior 
produces further consolidation beyond that already 
effected by the weight of earth material 

When the surface load is applied, the pressure pe! 
square foot at the upper boundary of the compressibl 
laver is 1000+ 3,240—4,240 pounds or 2.12 tons. That 
at the lower boundary is 3,460 or 1.73 tons. From 
figure 2, these pressures correspond to voids ratios o! 
0.90 and 0.94, respectively, the average value of e being 
0.92=e. Prior loading, the average value of e, as 
found in example !,:is 1.17. The ultimate settlement 
is therefore 

1.17—0.92 


a] x10 


r= 13.8 inches. 
. 1+-1.17 


1.15 feet 
The average \ olds ratio is reduced by the surface loading 
from 1.17 to 0.92. 
corresponding value of ¢ from figure 3 is 0.0138. 

Exrample 2 (continued).— Boundary condition A; tv 
stratum is sandwiched between sand layers. For the 
time-consolidation computations, again apply equatio! 
(17) for boundary condition A (see fig. 8) 


WAG) lasek) 


For g,;=25 percent, t — . 0.0.39 
1400 ¢ 1,400*0.0138 
year=12 days. 
Similarly, #=0.135 year 
0.325 year 


49 days for g,=50 percent ed 
119 days for g,=75 percent. 


These values are nearly equal to those computed ‘0! 
boundary condition A, example 1, and would be he 


The average of these is 1.045 and thie 





a tase 








ol 


he 


ee ede) ler Mesh Sie 





March 1937 


Y 000 1000 TYTN POITIER STITT ITI 


AD) MPER\ ‘ . +4 A frtity 


priate 


Figure 9.—-DistRIBUTION OF MaxiImMuM PRESSURE UNDER 


RECTANGULAR FoorinGc. Sor LAYER OVERLAID By SAND 
snp RESTING ON Rock. 
same were it not for small variations in ¢ as the pres- 
sure 1S varied. Kor most practical purposes, when 


boundary condition A obtains, the time-consolidation 
relation can be considered independent of the amount of 
pressure and the type of pressure distribution 

Example 2 (continued)-— Boundary condition B; the 
compressible stratum rests on rock and is overlaid with 
sand. When condition B obtains in this example (see 
fig. 9), use equations (24) and (26), together with tables 
2 and 3. The total settlement, Q, is the same for both 


boundary conditions, A and B. The value for uw (con- 
sider only the active pressure, fig. 9) is 2.22. Then 
J’ is obtained from table 3 and is taken as 0.67. From 
table 2, when g;=q3=— q;=25 percent, V,;—0.12 and 
V,— 0.03, and from equation (26), 

V,=N34+-J’(.N,— Ns) = 0.03 +- 0.67 (0.12 — 0.03) = 0.09. 
This value for .V; is substituted in equation (24) and 


tuking e=0.0138, then for q; percent, 


ND,’ 0.09 (4.61 9 909 36 d 
099 vear=st uys 
1400 ¢ 1400™0.0138 yn dls 
Koy a= s—q 50 percent, N, (table 2 0.49 and 
N.=0.24. 
Then N;= 0.24 +-0.67 (0.49—0.24)=0.41 and 
0.41 4.61 
0.451 vear=—165 days 
1.400 0.0138 tol yeal ‘allie 
™ ilarly for 5 » percent, 
1.08 (4.61)° 
9 yvears=435 davs 
1.400 0.0138 119 ¥¢ oe 
ample 3. Case 4, pressure distribution; type E 
lox ding; boundary condition B. It would seem that 
ciso 4 pressure distribution is a comparatively rare 


oc irrence. Problems of this type may be encountered, 
ho ever, when consolidation produced by the weight of 


eal material alone is incomplete at the time of load 
ap ication at the ground surface. Figure 10 illustrates 
suc a condition. If it is assumed that consolidation in 
th ompressible stratum 5 feet thick resulting from the 
We it of sand above is only 50 percent comple ted, then 
thi. \eight exerts an active pressure which is zero at the 


to al 1,000 pounds per square foot at the bottom of the 
coi_pressible layer. The distribution of active pressure 
res 'ting from the weight of sand alone is, therefore, 
triingular, It is assumed that consolidation of the 
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compressible layer as produced by its own weight is 
complete at the time of loading. 
Referring to figure 10, a load of 3 tons per square foot 
s applied over a rigid circular footing 10 feet in diameter. 
I is desired to determine the total settlement and the 
rate of settlement of the compressible layer which rests 
on rock and 1s covered with the sand layer 10 feet thick, 
50 percent of the consolidation producible by the weight 
of sand being completed at the time of load application. 
To compute the pressure, p,, at a point on the vertical 


centerline 10 feet below the rigid footing, substitute 
values in equation (12). The cosine of angle B (fig. 10) 
10 2 . rn 
Is = = and Pp equals 3 tons per square foot. Then 
yl2o yo 
) = 
p(i—cos’ps 3( 1 - =) 0.854 ton per square 
) 


54/! 
loot 

At a point 15 feet below the footing and on the verti- 
cal center axis, p,==0.438 ton per square foot. There- 
fore, the applied load produces a pressure of 1,707 
pounds per square foot at the upper boundary and 876 
pounds per square foot at the lower boundary of the 
compressible layer 

Prior to loading, ¢ at the upper surface of the layer is 

28, corresponding to a load of 1,000 pounds per square 
foot, and at the lower boundary it is 1.44, corresponding 
to a load of 500 pounds per square foot. The average 
is 1.36=e,. The voids ratios (see fig. 2) corresponding 
to pressures of 1,000+-1,707 pounds per square foot at 
the top and 500+-1,876 pounds per square foot at the 
bottom of the layer are found (fig. 2) to be 0.99 and 1.01, 
respectively. The average is 1.00=e,. The total set- 
tlement is, therefore, 


0 i 367 08 ck 
. l 1.36 


[t is important to note that prior to loading, soil at the 
top of the compressible stratum has reached equilibrium 
under the weight of sand (1,000 pounds per square foot) 
and that soil at the bottom of the stratum has not as yet 
“i en consolidated by the sand but has been ¢ consolidated 
by the we ight of compre ‘ssible soil, i e., 500 pounds per 
square foot. Therefore, the voids re tio at the top of the 
layer prior to loading corresponds to a pressure of 1,000 
and that at the bottom to a pressure of 500 pounds per 
square foot. After a very long period of time subse- 
quent to loading, the soil at the top of the compressible 
layer will be in equilibrium with a pressure of 2,707 
pounds per square foot and that at the bottom with a 


0.76 foot. 
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Figure 11.—-DistriputTion or Maximum PRrReEssSURE UNDER 
Strive OF INFINITE LENGTH Sou. Layer OVERLAID BY 
SAND AND RESTING ON Rock. 


pressure of 2,376 pounds per square foot. These pres- 
sures correspond to definite voids ratios, found experi- 
mentally (see fig. 2), and are used in computing Q. 
When load is applied, the active pressure at the top 
of the layer is 1,707 and that at the bottom is 1,876 
pounds per square foot. Hence u wide 0.91, a 
1.S76 
value between 0 and 1 (case 4). Apply equations (23 
and (25) and use the values for .V, and N, corresponding 
to different values for g;=q., table 2. Values of J are 
found from table 4, and for uvw=0.91, J=0.06. For 
(i= qz=25 percent, N,;=0.12, N.=0.31, and 


N,=N,4+-J(N.—N,) =0.12 + 0.06(0.31—0.12) =0.13 
] - ‘os 
D, K 9 =sagn X9=2.12 feet. 
1+ é, 2.00 
For e= ~~ —=1.18, c=0.0106 (fig. 3 


, N,De 0.13 X (2.12)? 
Then for g,=25 percent, ¢ Nee enn a 
0.039 vear= 14 days. 
ev ges - 0.50 * (2.12)- 
Similarly, for g,=—50 percent, ¢ 1.400 <0.0106 O.15 
vear=55 days, 
1.19 X< (2.12) 


and for qy=75 percent, ¢ 0.36 year 


1,400 « 0.0106 
132 days. In 132 days the settlement is 75 percent 
of 0.76 foot or 6.8 inches. 

Example 4.—Case 5 pressure distribution; type D 
loading; boundary condition B. All other conditions 
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ure the same as in example 3 except that in this case 
the rigid footing is a loaded strip 10 feet wide and of 
relatively great length instead of a circular bearing area 
and consolidation caused by weight of sand is complete 
prior to loading. Figure 11 shows a vertical cross 
section. To find p, at the top of the compressible 
layer and at a point on the vertical center axis through 


the cross section, apply equation (11). The angle, a, 
is expressed in radians and sin a=0.80. 
Pp , 6,000... 
P: a+simMe ‘ 14 (0-98 0.80) =—3,290 pounds per 
Tv » . 


square foot 
At the lower boundary and on the centerline, 

6,000 ; — ; 
p 2 14 (0.64 + 0.60 Zucto pounds per square foot. 
This is case 5 pressure distribution, and u—1.38. The 
average voids ratio prior to loading is the average of « 
corresponding to 1,000 and of ¢ corresponding to 1,500 
pounds per square foot and is 1.2l1—e,. Subsequent to 
loading the values of e at the top and bottom of the 
compressible stratum correspond to the equilibrium 
pressures, 4,290 and 3,875 pounds per square foot, 
respectively. From figure 2 these values of e are 0.90 
and 0.92, the average being 0.91—e,. The total settle- 
ment is therefore 


1.21—0.9] F 5 ; 
4 | 1+1.9] ] <5—0.68 foot or 8.2 inches 


To compute the rate of consolidation, apply equations 


24) and (26) and make use of tables 2 and 3. J’=0.87 


L.ObL¢ 


» 


for u=1.38 (by interpolation). Fore 


0.0134 
- N;Dy? 0.11 X (2.26) 
For (1 =4s=%s=25 percent, f 1,400e 1,400 0.0134 
0.03 year=11 days. 

0.46 * (2.26)° 


1,400 0.0134 0.125 vear—46 


For g;=50 percent, f¢ 


days. 
L.14 X< (2.26)° 


O310 ar ” 
1,400 « 0.0134 l veal 1] 


kor 5 49 percent, t 


days. 
In 113 days the strip settles 75 percent of 8.2 inches 
or 6.2 inches. 

The preceding examples are intended to illustrate the 
application of the theory of soil consolidation. Part 2 
contains the derivation of equations and the develop- 
ment of the basic theory. 








eee 
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PART 2.—DEVELOPMENT OF THEORY OF SOIL CONSOLIDATION 


In 1885, J. Boussinesq ® developed in considerable 
detail and usable form a solution of the general problem 
of determining the distribution of stresses and strains 
in a semi-infinite, elastic, isotropic solid, bounded by a 
plane and loaded by a single concentrated force at a 
point on that plane. The formulas of Boussinesq have 
had wide application in the development of the science 
of soil mechanics. 

The distribution of stresses across any plane parallel 
to the plane surface is shown in figure 12. If a single 
concentrated vertical load, P?, is applied at a point on 
the surface, then at a radial distance, R, from this 
point and at a depth, z, below the surface, the stresses 


ut the point (r, y) are: 
Dr an R parallel to PR 30 
p oo Re parallel to v axis 333 
Py ann = parallel to y axis (34 
3P 


parallel to 2 axis 


I: Dah? R* 


The derivation of these equations may be found in 
almost any standard text dealing with the theory of 
elasticity.” These equations contain no elastic con- 
stunts, however, and for this reason it has probably 
been assumed that the distribution of stresses does not 
depend upon the type of material. Cummings,‘ 
Krynine,® and others have discussed this point at length. 

A complete analysis of the stress distribution includ- 
ing planes not parallel to the boundary involves 6 com- 
ponents of stress, 3 shears and 3 normal stresses, and 
the elastic constants appear in some of these equations. 
In the theory of consolidation of soils we are concerned 
tainly with the vertical normal stress, p,. At a given 
depth, z, it is seen from equation (35) that p, is greatest 
when R is smallest, that is, the maximum values for p, 
ure along the vertical centerline of the loaded area 

x=y=0, fig. 12). 

The value of p, at any 

wing that R?=r?+2?=r*+y?+2?. If a combi- 
nation of point loads is given, we may obtain the stresses 
ut a point resulting from each load considered sepa- 
raicly, then add these and obtain the total pressure at 
the point produced by the combination of loads 
(uiboy '® has shown that by dividing a uniformly 
ded area into rectangular elements so that the 
I 


point is easily computed, 


onger side of any element is less than half the distance 
rn the center of the area to the point, the computed 
Valle for p, at the point will not deviate from the 
co-tect value by more than 6 percent. The accuracy 
is -reater with further subdivision of the rectangular 
e ents. 

lication des Potentiels a l’Etude de I’ Equilibre et du Movement des Solides 
1eS. Paris, 1885 

example, see chapter 1! of Theory of Elasticity, by S. Timoshenko 

tribution of Stresses under a Foundation, by A. E. Cummings, Proceedings, 
in Society of Civil Engineers, vol. 61, no. 6, August 1935 ' 
cussion of Cummings’ paper by Clement C. Williams, D. P. Krynine, and 
Wileoxen, Proceedings, American Society of Civil Engineers, vol. 68, no. 8, 
pt October 1935 

irths and Foundations. Progress Report of Special Committee, Proceedings, 
can Soctety of Civil Engineers, vol. 59, no. 5, May 1933 
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FIGURE 12 DISTRIBUTION OF STRESSES ACROSS ANY PLANE 


PARALLEL TO THE PLANE SURFACE. 


FORMULAS DEVELOPED FOR DETERMINING PRESSURE AT POINT 
BELOW LOAD 


With reference to figure 12, if the expression for 


P2=>5s—p; JS iMtegrated in the y direction from 
' 2rh° : 
y © to y=+o and the resulting “line load”’’ is 


multiplied by dz, we obtain the total load on a strip 
infinitely long and of width, dz, produced by the point 
load, P, at the boundary. In this integration, z and z 
are constant. Thus: 


3Pz 3P2 


Is e _ ? > s 6360 RID 
Bs 2nrR 2a(z+y +c) 


‘ 3P2 (°T 2?4+7y'+ 27 75, a 
p da J ies *dy 
T 0 ieee pon . 
This is easily integrated if we put y Via tan 6 
see fig. 12). Then 
q 3Pz were 
Jp. (da = = a2 | cos*6dé, 
or 
. 3P2 sin®6 “]*/? 2P2 
. = <3] Sin @ a dx = ae 36 
| Pe (ae >a 3 I raat? W--(36) 


If this expression is again integrated in the x direc- 
tion, taking twice the value obtained between the limits, 
x=0 and r= & (fig. 12), the result is equal to P, the 
applied force. Obviously, on any plane of infinite 
extent and at a finite distance below and parallel to 
the loading surface, the total distributed load is equal 
to the applied load at the surface. The mathematical 
result is in agreement with this self-evident fact. 
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In equation (36), if P is taken as p, the vertical pres- 
sure per unit area, and if p is considered as concentrated 
at the center of a unit area, then the total vertical 
pressure on the strip at a distance below the surface is 
2pz 

dx. 


equal to a This is the total!pressure, /‘p., 
w(2*-+ 2°)* - 


on a strip below the surface produced by a unit load on 
the surface. Conversely, the total load on an infinitely 
long strip at the surface and parallel to the y axis, 
produced by applying the unit load, p, over the length 
2pz 


of this strip will cause a vertical pressure vo da, at 


Ts K 

a point below the surface. Briefly, a point load at the 
2pz 
w(z°+ 2°) 
below the surface and conversely, a strip load of unit 
load per unit area on the surface produces a pressure, 

Qn 
ae 5dr at a point below the surface, action and 
a(z°-+ 2°) 
reaction being equal. It will be noted that the integra- 
tions of the Boussinesq formulas are with reference to 
areas or points below the loaded area at the surface. 
At the surface, z=0 and there is nothing to integrate 
that has reference to the fundamental equation (35). 

Figure 13 is a vertical cross section of a uniformly 
loaded area (p per unit of area) of infinite length and 
of width, 2b. The load on a strip of infinitesimal width, 
dz, produces a pressure, p,, at a point T equal to 

2pz2° 

—* dr. 
We 2)° 
limits, there is obtained the value for p, at T resulting 
from the total load over the whole area, infinite in 
length and of width, 26. The integration is more easily 


surface produces a total load, 5dr on a strip 


If this expression is integrated between 


*an J . ry 
performed if we let _=cot @ (fig. 13) Then at the 


point T, 


» 2 2\ -2 ») "ay 
2b f* (+5) de? 
Te z-—<2b 3 T Ja 


) ‘ 
=Pig -SIN @ COS(a, + a,)] (3 
T 


(1 + cot) ~“ese*adé 


Pz 


“J 


Y 


If the point, T, is on the vertical centerline through 


the vertical cross section (fig. 13), then a»=180°—a, 
and equation (37) reduces to 
pP . 
Pz (a-+sina) (38) 
us 


For uniform pressure distribution (p per unit area) 
over a circular bearing area of radius r (see fig. 14) the 
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pressure at a point below the loaded area and on the 
vertical centerline is, 
"ts ("2* 3p : . “ rdy 
P | | Jr in rdrdé a y (7 
Integrating and simplifying, 
| ; , 
Pz Pp l ; and since tanp 
( et ) 
| 
p pa cos’ p OY 


MECHANICS OF SOIL CONSOLIDATION OUTLINED 


Consider the flow or egress of water from a wet, 
compressible, soil stratum subjected to vertical loading 
The soil is relatively impervious and is ‘‘sandwiched”’ 
between two permeable sand-gravel layers, L, L, figure 
15. The movement of water saturating the stratum 
is in both the upward and downward directions, 
parallel to the vertical axis, ZZ, and directed toward 
the two permeable layers, L, L. 

At any finite time, ¢, after load application, that part 
of the total pressure carried by the water, p,., is consid- 
ered as the ‘hydrodynamic excess.’’'' The remainder 
of the load is carried by the solid material, the “soil 
skeleton.”? At the outset of loading (f=0) all of the 
pressure is carried by water if the compressible layer 
is saturated. At any subsequent time, the hydrody- 
namic excess, Py», is zero at the surfaces of contact o! 
the soil stratum and the permeable layers, but elsewhere 
it is a finite quantity, varying with the depth, z. The 
upper and lower faces, AB and CD of the elementary 
prism, ABCD, are each of unit area. The thickness 
of the elementary prism is Az (fig. 15). 
| We may then define the following expressions: 

p=the total pressure per unit cross-sectional ares 
lt includes the pressure on the soil skeleton and the 
pressure, P,, carried by water alone. 

Pw=the pressure per unit cross-sectional area that is 
carried by water alone. 

P— Pwo= the pressure per unit cross-sectional area that 
is carried only by the soil particles, i. e., by the soil 
skeleton. 

a=the rate of change of the voids ratio with respect 
to (p—p,) and is called the coefficient of compressibility 








ko=the reduced coefficient of permeability; see equa- 
tion (4), Part 1. 

dp, : ; : 

: =the rate of change of the hydraulic pressu’e, 


Pw, With respect to the depth, 2. It is the hydraulic 
pressure gradient. 


Lp, Lp, 
tf 


. » C ° ! 
2 =the rate of change of / with respect to 2 and 
zZ az 


1! Terzaghi, K., Erdbaumechanik 


Franz Deuticke, Vienna, 1925 
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is taken as constant through the very small vertical 
distance, Az 


dp. _ dp ; ‘ 
az—F"—the total change of “2” within the inter- 
dz d: 
val, Az. 
dp . 5 
ke j- Az—the rate at which the elementarv prism 
det 
loses water. It represents the difference between the 


quantity of water flowing from the upper face, AB, in 
unit time and the quantity flowing into the prism 


through the face, CD, in unit time. It is therefore 


the rate of volume reduction of the elementary prism 

dp d p Ds . . . 

- the rate of change of p—p, witl 

dt dt , 2 ee ee 
respect to time 

dp, ; 

a it the rate of change in voids ratio with respect 
to time 

dp 2 . 

a it Az=the rate of change in volume of the ele- 


inentary prism with respect to time 
dp, dp, oP 
P me Fl dt 
the partial derivatives. The faces, AB and CD, of the 
elementary prism are both of unit horizontal cross- 
sectional area and are parallel to each other. 

Che amount of water lost by the elementary prism 
a 
" a2" 
tween the amount of water flowing out from the face 
AB and thet flowing into. the prism through the face 
Cl). But the rate of loss of water (which is the rate 
of reduction in volume) by the elementary prism is also 


The expressions are considered as 


in unit time is hf Az, which is the difference be- 


dp, ee 1? p, dp —— 
equal to a 5 Az Chen koa aZz=—d F A Divide 
bot sides of this expression by a and write ¢ for 
Then, sinee Az cancels, 
LP dp, 1() 
az dt 


| ‘ils is analogous to the development of the equation 
lor ‘he flow of heat. A more complete derivation otf 
eq! ion (40) has been published.” 


rN CONSOLIDATION RELATIONSHIPS FOR UNIFORM 
DISTRIBUTION DERIVED 


PRESSURI 


olution of equation (40) is 


Pw=e-k*(A cos Kz+B sin Kz) 41 
. wl 


e A, B, 


: thematical Theory of the Process of Consolidation of Mud Deposits, by 
Bee Ortenblad. Journal of Mathematics and Physics, vol. 9, no. 2 


and AV are constants. 
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In figure 16, the reduced thickness (the thickness the 
compressible stratum would have were there no voids 
between the particles) of the soil laver is 2/) where J, is 
the maximum distance that water travels if the layer 
That is, dy =2l. Water saturates 
\pply load. At the very instant that load 
is applied, as noted before, the water in the soil takes 
all the pressure (see fir. 17—A That is, at ¢ =), the 
total pressure, p, is equal to p, at any depth, z (see fig. 
I7-A At any other time after loading, p, is less than 
the total pressure, some of which has been transferred to 
the soil skeleton. The condition, subsequent to apply- 
ing load, is represented by figure 17—B. Atany level, z,a 
part of the load, p,,, designated as the “hydrodynamic 
excess,’ is carried by water. The load carried by the 
soil skeleton at depth, 2, is the total pressure minus p,,. 

\t any time subsequent to f=0, p, is equal to 0 when 

0. The water escapes instantaneously at z=0. 

The same condition prevails at the lower boundary, 
that is at 2l), Po=0 if t is greater than 0. At t=0, 
Pe is a function of 


has a thickness, d 


the soil 


For present purposes, assume that the total pressure, 
p, (that on solid and that on water taken together) is 
uniform and constant throughout the depth, z. 

We now consider equation 


41) under the following 
boundary conditions 


Bsin Kz)e x 
Condition 1: p,—0 when 0 and t>0. 
The sine of 0 is 0 and sin Kz=—0. 
Therefore, since p,—0 and B sin Kz=0, it follows 

that A cos Ke must also =0. But cos 0—1 and 

Acos 0=A. Then A must equal zero. So condition 1 

reduces equation (41) to 


Pe (Acos Kz (41) 


p Bsin Kz) e-ck*s (42) 
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Condition 2: When z=2/, and p,=0, then sin Az 


sin 2K, = Q, 
Hence Kz must be an integral multiple of 7. 


z=2l, and Kz=2Kl, for condition 2, and sin Kz—0) 


when Kz=2Kl=7, 27,327...nn, where n is any 
integer from 1 to infinity. 
9 ? 
+e > » us at T ais 
“oD ° 

Since 2Al,=nz, K of, 91, 91, = 91 

mn ° nénct : P 

The expression, e ‘**, becomes e¢~ 4: Substitute 


Qn nT 


for A the suecessive values = ‘te. in equati 
or A the successive vaiues, 2], 5], 9] yete. in equation 
(42), add and obtain: 


réct 


4r2ct 


; - a5 - FE dle tn 2 FA 
Po= Ue sins + ast sin) + 
— 3] am Ot) 
néwéet } Es 
" tT 
+a,e “” sin= (43) 
2l, 
When t=0, equation (43) becomes: 
«88 . 282 _s 372 
) Pu=p=Q sing +42 sin af, 1% Sin Gy + 
» REZ 
; +a, sin 9] (44) 
r This is true for all values of z between 0 and 2/,. To 
i 


. . « FZ ; 
find a,, multiply through by sin5y and integrate all 
alo 


terms between the limits, 0 and 2/,. Then, as usually 
happens in the case of a Fourier series, all terms on the 
right, except the first, vanish, and since p,,—a constant 
=p when t=0, 


SYS? 2 SSeS SoS 


*2lo ~ 2 i 
- We ° T 
then | sin=>dz—=a [ sin?—dz__ 45) 
P 0 2l * 0 2l, a 
i 4 4p 4p 
( ay = ’. (l9= 0, a,—= - I ;é¢.= / ’ 
Sn nar 
} where n is odd. 
Substituting the values of a,, a2, a3, -* +a, in equa- 
: tion (43), there results 
i 
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4p rict rz, 4p Oxict a 
rh oe di? oo ,U** 
| ) e sin rr ~ sim pattie 
| Pw 7 D1, 3m 21, t 
2n+1)2e%ct 
4p was? . (2n az 
‘ 4lo? 46) 
(2n+1 im mai 21, a6 
Where n is 0, 1, 2, 3, ete. 


Multiply equation (46) through by dz and integrate 


between the limits, z-=0 and z=2/). We have: 
(nade o it [inEZaeet we nrg 
2 o = Side +--e ™ SINS) dz 
0 ! ’ T J0 21, oT J0 21, 
2n+1)*rict po 
4p — ¥ - (2n+1)xr2 
410? . ed - gf 
+ a? e sin d: (47) 
(27 I l)r 0 21, 
*2/ Slop rel r: 2/ 
or ) dz=— e@ 4? | COS 
|, Pu 2 2ly ly 


Spi, we nz | 
re lo? | COS . 
ar* 21, 
Spl wash) ann (2+ Lise i> 
= 1 Alo? COS : 
2n-+ 1)*x* 2, 0 


,D ‘ ‘ “n+ 1)*s 
16) LN l és 7 
Tr Pe (2n-+-1)>- 


n 0 


This ts the total p, at any time greater than zero. 


The total p, at ¢ 


) | Dp d: Pleo? 2lap 
since Py=—p when to. 


. 


p »dz—total p, at any subsequent time, f 


So" Pre Az 
2lyp 
time, ft. 


the percentage of consolidation at 


15 


That is 


16 pl, : e @-a¥ «6 
ae ) 
x )) 33 49 
7 2 pl 
/ ON »—25N y FON 
’ —s r(¢ 4-4" + + -) (48 
or q; 9 95 49 
, wet 
where N=-—,, 
41, 
Figure 18 illustrates diagrammatically the load 


fractions carried by the clay skeleton and the water. 
In figure 19, X = load fraction carried by water, and 
Y load fraction carried by soil. 


The total p,= |" Po dzand X+Y 


ae dz 


total pressure. 


i Pw dz 





y 
-]— “5 ==> when t=?, >0 
1 joy ae os hen 
¥ wict 9rict ] (2n+1)*r'ct 
0;=1 —8, r( c ~ 4p? Ge ~ 413 see — Ab) “4 
fi t + (2n+1)* 
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FicurRE 19.—-TRiIANGULAR PRESSURE DISTRIBUTION. THE 
LOWER TRIANGLE IS THE MIRROR IMAGE OF THE UPPER ONE. 


DERIVATION OF TIME-CONSOLIDATION RELATIONSHIP FOR TRI- 
ANGULAR PRESSURE DISTRIBUTION, PRESSURE INCREASING 
WITH DEPTH 


lor a triangular pressure distribution (see fig. 19 
the two boundary conditions, 


Condition 1: p,=0 at z=0 


9. IPw 


and condition 2: :(0 at z=/ 


ist satisfy equation (41), 


Peek (A cos Kz+B sin Kz 


from 


i) 


condition this equation reduces to equation 


as before and 
Po= ek (B sin Kz) --- — (42 


r the second condition, at z=, 


dp, . 
ai : =e-k*« BK cos Kz 
d- 


O0—e-K*« BK cos Aly. 
mar 
= 


cos Al,=0 when Kl, 


1 is any odd integer. That is, 


- w on or mr 
K 


eee rar ace 6.69 Sane 
2ly 2ly 2, 2ly 


Ss 


itute these values of A successively in (42) and 


adc the results together; we thus obtain 
wet 9rict 


~ 


ae P — - ‘“ 
: , » 4le o as 4le? ¢ 
} Pw=aye ’ ma tae * sin 


ww 


us 


a 


to infinity. (49 


N 


At , equation (49) becomes 


312 - Mrz 


I +--+ +d, sin [i (90) 


° T2 . 
a; SIN 5> +42 sin 
2l, 











mm Whe is 1, 3, 5, ete., i. e., any odd, positive integer. 
bi | oeflicient, dp, 1S ev aluated in the usual manner. 
| At t=0, py=f(z)= n between the limits <=0 and 
0 
m 2=l. At t=0, between the limits z=, and 2l), 
“on 
Pw | 2 on as 
ial 





0 
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FiGurRE 20.—TypEes OF PRESSURE DISTRIBUTION IN A Com- 


PRESSIBLE LAYER. 


| bon z 
Um il f ‘, sin si = de Te dz+ | (a? 


MZ 
re) si n> = dz | 
“lo 
Where m is any odd, positive integer. 


Sp, e —* 
PF. for m=), 5, 9, 13, 17, etc., 


Then 





Integrating, @», 


mT 
and 
Sp » — = 
n= ——, for m=3, 7, 11, 15, 19, ete. 
he 


By substitution equation (49) becomes 


ret QOxct 





SP; 4l2 . WZ Spi Aly? 37Z i 
, sin ; e sim—=-it... ----- (51) 
rT” > Ya 21, 
m m—1l 4 mix? 2ct lo 
Au 8n,\ \(—1) 2 4lo? . Mrz 
Then | p,.dz= An) ——— e sin 7 dz 
—«. m? 2, 
ea 0 m= | 0 
ict Oxrrct . ‘~. (2n+1)*x%ct 
l6pyl, ae? _16pilo rio 16 plo” ._ Go? 
~ ° 27% - *** (Qn-+1)3x* 
Where n=0, 1, 2, 3, ete. 
Piz “le rl 
\t t=0,f and} f lz PY 
l, J0 “ 
* lo 
Pd 
. 0 
( war 
f(z)dz 
JV 
J t 1 Ox t 1 25x ct 
oZ = - 
l aaa e 4 2 : 4 fi ’ . 4lo 
r Th + i595 + 
af -? 
rct 


32 » § ee ] a ‘ 
@=—1— 2(e- ——e % +755 Sete oe .) --(62) 


FORMULAS APPLICABLE TO OTHER TYPES OF PRESSURE 
DISTRIBUTION DEVELOPED 


General relations.—Figure 20 illustrates five different 
types of pressure distribution. In this figure the con- 
ditions illustrated are for a compressible stratum over- 
laid with sand and resting on impervious rock. 

The percentages of consolidation for the five types of 
pressure distribution shown in figure 21 are expressed 
as follows: 
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FIGURE 21.—VALUES OF g AND N For DIFFERENT VALUES OF u. 


Case 1.—Uniform or rectangular pressure distribution. 
S ae Wee l re 
n=1—(e*+e¢ Wor...) 

T \ ‘ <0 / 
Case 2.—Triangular pressure distribution, pressure = 0 


at upper surface. 
sea f l ; l 
- N —9N N 
o= —_ — e — =—@ eV, +...) 
d2 =l\ 7° T1235 


Case 3.—Triangular pressure distribution, pressure = 0 
at lower surface. 


o(a—%) 


q3 D =2q1— 2 
2 


_ Case 4.—Trapezoidal pressure distribution, pressure 
increasing with depth. 


DP 

PH T 2 9 

“.=- ol 
p+P 


Case 5.—Trapezoidal pressure distribution, pressure 
decreasing with depth. 





) 
py— q2 (p—pi)gt* 3 
qs>= —" De —-= 
) 
’ ai “ p—— 


If the compressible layer is in contact with permeable 
sand layers above and below, then equation (48) holds 
for all of the five types of pressure distribution shown 
in figure 20 and interpolation formulas are not needed 


PUBLIC ROADS \ 





1. 18, No. 1 


in this comparatively simple case for, if q, the per- 
centage consolidation, is plotted against the time, the 
curves for the different pressure distributions will 
theoretically coincide for the same soil. The ultimate 
amount of consolidation may vary with the type of 
pressure distribution but the percentage of the ultimate 
consolidation attained at the end of a definite period 
of time is theoretically independent of the pressure or 
its manner of distribution if the relatively impervious 
stratum is confined between two relatively permeable 
lavers, and the rate at which gq; increases depends 
solely on the permeability, compressibility, and thick- 
ness of the clay layer. 

The expressions for cases 4 and 5, trapezoidal pres- 
sure distributions, are combinations of equations (48 
and (52) and are obtained by integration according to 


the procedure followed in deriving equation (52 The 
expressions are, case - 
dy 
PUT 5 8 
4 - ays) 
,?P 
PT 5 
and ease 5, 
PY =e 
q ve D4 
Py 
P~>5 
Case 3, triangular pressure distribution, heiny 


1 
infinite, is a special case of equation (54 For case 


pp, and hence from equation (54), 


q3= 2M — 42 de 
When the expressions on the right of equations (4s 
and (52) are substituted for 1 and q2, equation a 
is written 
16 To , i 
q 2 37 8” ae si T Sree : 
T \y yxy) _j 
320 it , 
1- e~* —@—9N =e ; a6 
24 125 
For mulas for tra pe ~oidal pre SSIUTE¢ distr ihut iO? i. Witl 


reference to figure 20, both q, and q; approach q as : 
limit as p,; approaches the limit zero. This is also 
evident from equations (53) and (54). Note also that 


‘ Die . . 
p+ is the area of the trapezoid, case 4, when / 1 and 


1 


) . . . ~ . 
. o is the area of the trapezoid, case 5, for / I 


Pp 
For the condition, ¢q.=q=q, the numerical valu 
N,, will be between that for N; and .N2, that is, 


N,=N,+Jd(N2—M™)). ea 
ex. , — - Es 
This is equation (25), Part 1. Then ie 

(1—q,) =(1—q.) = (1—q,) and log, (l—¢ ee | 
log, (l1—q2) =log, (1—q,) =log, R. 
Neglecting all terms in e except the first in equations 
(48), (52), and (53), 
Ss , --\ 
log, (l—q) =log. =—MNi=log, R_-_----9' 
1 ae 
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32 y ~ 
log, (l—q.) = log, 3 —N2 log, R - (58) 
Ss 6 
rl r vsP: 
log, (l—q,) =log, N,=log, R __(59 
- } 
pty 


where log, denotes the natural logarithm. Equations 
(58), and (59) are easily solved for .V,, N2, and 
N, and the values are substituted in the expression 


(od J, 


] 


- N,—N, 
A = - (60 
N, N, 
which is equation (25) solved for J. Thus 
») } 
log | 2 aps | 
J wilp - P1) 61 
log = 
T 
, u ; 
Since for case 4, u P or p ms » then by substi- 
PP l Ue . 


tution in equation (61) one obtains 


| 2 2)u- 
or 

w(l-+-u 
J 


, 62 
log 
which 1s equation (21), Part 1. 
For the condition, g;=qs=gs, one May write 
y’ N;—Ns3 2 
. = ; a >) 
N,- N; 


Equation (63) 1s equation (26), Part 1, solved for J’ 
furthermore, 1—q,;=1—q;=1—2q+q, from equation 
55), =1—qs. 


By dropping all terms in e except the first, 


8 . ; sie 
log, (1 log,5— N,=log, R as before. __(57 
a 


-qi) 


log, (]- 3) 


loz 8-95 |- Nz log, R 
us v 


—q; may be substituted for gq: 


64 


equation 
to obtain 


94), 24, 


) 
(p—p)ats a 


q ees | 
p—f 
ox 8 ks ( 8 4 
2s as 
log. (1- q;)=log,| - = re N; (66 
red 


Alter substituting the values obtained by solving 
equations (57), (64), and (66) for N,, N3, and N; in 
Pr 


" me 2 \ ‘ ‘ U 
equation (63) and substituting - 7 for p, the expres- 


sion obtained for J’ is 
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u(rm—2)+2 
log ( 2) (u+1) 
r—?)(u+ 
J : 


- 67) 

T 

log] - 
2(x—2) 
This is equation (22) of Part 1. 
1 

The expression for J was derived by Gilboy.” The 
expression for J’ was given by Kimball." Neither 


Gilboy nor Kimball indicated the derivations and both 
have given negative signs for equations (62) and (67). 
It would of course be impossible to obtain table 3 of 
Part 1 of this paper with negative values for J and J’. 
Since only the first terms in e were used in their deriva- 
tions, the expressions for J and J’ are of course only 
approximate 

Figure 21 shows curves obtained by plotting values 
for N on log seale as abscissae against corresponding 
values for q plotted as ordinates. The different curves 
are for 8 particular values assigned to u, these being 
u=0, %, %, 1, 2, 3, 10, and For intermediate values 
of u, values of g and N may be found by interpolation. 


CHARACTERISTICS OF THEORETICAL TIME-CONSOLIDATION CURVE 
DISCUSSED 
Equation (48) is one of a smooth curve having several 
important characteristics, one of which is illustrated in 
figure 22-A and 22-B. Figure 22—A is the curve 
; : ; ~ ¥'C . P 
obtained by taking t=—.N, that 1s, as unity In equa- 
tion (48 Let g, be the percentage consolidation at the 
end of a definite time interval, ¢, and ¢, the time re- 


quired to produce half this percentage of consolidation, 


L.e., = Qh 22-A, t=1.4 


For q,=—80 percent, figure 


c : ‘ . a 
minutes and the time, ¢;, required to effect = of 80 per- 


cent consolidation is seen to be 0.3 minute from the 
time consolidation curve, figure 22—A. The ratios, 


t ‘ 
figure 22-B. 


} 


are plotted against q,; in Within the 


a ; eee 

range q,—0 to q,=70 percent, the ratio, ~ is fairly con- 
stant and is approximately equal to 4. This indicates 
that within this range the time-consolidation curve is 


approximately a parabola. For values of q, greater 


a : fx ; 
than 70 percent, the ratios, ,? increase. An ordinate 


1 


value of q,=—90 percent on the theoretical curve, 


figure 22—A, corresponds to a value of -5.27, figure 


1 
22-B. The value, f=2.1 minutes, corresponds to the 
t ani 
value, $37. 
t; 
Since the relation between q, and _.N, is fixed, the ratio, 
t 


90 percent is invariant for all values 


5.27 for q, 


TC 
4l, 

As a rule the time-compression curves for relatively 
impermeable soils have very nearly the same shape as a 
theoretical curve up to a point which, on the theoretical 
curve, has an ordinate of 90 percent. Beyond this 
point the two curves diverge appreciably. Since an 
13 See equation (31) in ‘Earths and Foundations,” 


of Civil Engineers, vol. 59, no. 5, May 1933. : 
4 Proceedings, American Society of Civil Engineers, vol. 59, no. 6, August 1933. 
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| ordinate of 90 percent on the theoretical curve corre- 


=< ae ee a 
sponds to a definite value of Zo 9-27 , it is easy to com- 


1 
| pute the coefficient of consolidation of a soil by adapting 


the actual curve to the theoretical one according to the 


| following procedure: 


1. Plot the laboratory time-consolidation curve, 
corresponding to figure 22—A. 
a ; oe -_ 
2. From the time-consolidation curve plot the; ratios 
1 
corresponding to definite percentages of consolidation. 


. t . 
3. From the 74? % curve find the value of q, corre- 
1 


; t _. 
sponding to 75.27. 
1 
4. From the time-consolidation curve find the value 


; Ft ‘ 
of ¢ corresponding to the value of 4, for 7, = 5.27. Cali 
l 
this particular value ¢,. 
5. Substitute the value for ¢, thus found in the ex- 


; rc 
pression, 773 t,==2.1, and solve for c, 


214k? 0.851? | 
———— -=—= -—_—— ia ees eae. Bea (OS 
at. t. 


> 
. - 


c 


In testing with the Terzaghi compression device, 


the sample is confined between two permeable layers 


| 


] : 
and hence =35d. In this case, 


0.85d,? oe 
Ce gc oe cscecn nse =( OY 
This is another method of computing c. [See equa- 
tion (8).] 


ie 
a 
4 
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PUBLICATIONS of the BUREAU OF PUBLIC ROADS 








Any of the following publications may be purchased from 
the Superintendent of Documents, Government Printing Office, 
Washington, D.C. As his office is not connected with the 
Department and as the Department does not sell publications, 
please send no remittance to the United States Department of 
Agriculture. 


ANNUAL REPORTS 
Report of the Chief of 


5 cents. 


the Bureau of Public Roads, 1924, 


Report of the Chief of the Bureau of Public Roads, 1927. 


) cents. 
Roads, 


the Bureau of Public 1928. 


=. 


Report of the Chief o 


5 cents. 
} Report of the Chief of the Bureau of 
10 cents. 


. Report of 


10 cents. 


Public Roads, 1929. 


the Chief of the Bureau of Public Roads, 1931. 


. Report of the Chief of the Bureau of Public Roads, 1933. 


) cents. 


+ 


the Bureau of Public Roads, 1934. 


: Report of the Chief o 
» 10 cents. 
4 


; Report of the Chief of the Bureau of Public Roads, 1935. 


2 cents. 

4 Report of the Chief of 
10 cents. 

DEPARTMENT BULLETINS 

® No. 583D . . Reports on Experimental Convict Road Camp, 

} Fulton County, Ga. 25 cents. 


the Bureau of Public Roads, 1936. 


one ee | 





4 


® No. 1279D . . Rural Highway Mileage, Income, and Expendi- 


4 tures, 1921 and 1922. 15 cents. 
4 TEC!INICAL BULLETINS 
oy No. 2601 ... Electrical Equipment on Movable Bridges. 


35 cents. 


MISCELLANEOUS PUBLICATIONS 
me No. 7: 


“ 


Rock. 


25 cents. 


>) 


IP . . The Results of Physical Tests of Road-Building 


Federa! Legislation and Rules and Regulations Relating to 
Highway Construction. 15 cents. 
No. 191 .... Roadside Improvement. 10 cents. 
The Taxation of Motor Vehicles in 1932. 


An Economic and Statistical Analysis of Highway-Construction 
Expenditures. 15 cents. 


35 cents. 


Highway Bond Calculations. 10 cents. 








_ Single copies of the following publications may be obtained 
from the Bureau of Public Roads upon request. They cannot 
be purchased from the Superintendent of Documents. 


SEPARATE REPRINT FROM THE YEARBOOK 


No. 1036Y .. Road Work on Farm Outlets Needs Skill and 
Right Equipment. 


TRANSPORTATION SURVEY REPORTS 
Report of a Survey of Transportation on the State Highway 


System of Ohio (1927). 


Report of a Survey of Transportation on the State Highways 
of Vermont (1927). 


Report of a Survey of Transportation on the State Highways 
of New Hampshire (1927). 


Report of a Plan of Highway Improvement in the Regional 
Area of Cleveland, Ohio (1928). 


Report of a Survey of Transportation on the State Highways 
of Pennsylvania (1928). 


Report of a Survey of Traffic on the Federal-Aid Highway 
Systems of Eleven Western States (1930). 








A complete list of the publications of the Bureau of Public 
Roads, classified according to subject and including the more 
important articles in Pustic Roaps, may be obtained upon 
request addressed to the U. S. Bureau of Public Roads, Willard 
Building, Washington, D. C. 
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